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SUMMARY 


Equilibrium and frozen nozzle-flow properties for carbon dioxide are presented 
in the form of charts. In the calculations for the frozen flows the flow was assumed to 
be in equilibrium to specified points in the nozzle where the chemistry was assumed to 
freeze instantly. Then, the flows were further expanded with the energy in molecular 
vibrations either frozen or in equilibrium with translational temperatures. Results 
are presented for total pressures ranging to 1000 atmospheres and total enthalpies 
ranging to 55 kilojoules per gram. The properties of temperature, pressure, density, 
velocity, dynamic pressure, Mach number, Reynolds number, molecular weight ratio, 
and species concentrations are presented in charts. Temperature, pressure ratio, , 
and density ratio across normal shock waves in a nozzle are also included. j 


INTRODUCTION 


All the suggested model atmospheres for Mars and Venus contain varying 
amounts of carbon dioxide. In the absence of definitive estimates of gas composition, 
it is appropriate to investigate flight in carbon dioxide as one limit of possible atmos- 
pheres. To do this, it is necessary to document its thermodynamic and flow proper- 
ties. These properties for carbon dioxide in molecular form, given analytically in 
reference 1, correspond to low flight speeds. For total enthalpy levels in carbon diox- 
ide corresponding to higher flight speeds, the equilibrium thermodynamic and one- 
dimensional flow properties have been presented in graphical form in references 2 
and 3, respectively. In addition, the properties of flows across plane shock waves in 
carbon dioxide for speeds to 16 km/s have been presented in reference 4 for assumed 
equilibrium flows and for flows with frozen chemistry and molecular vibrations either 
not excited or fully excited. 

Not all nozzle flows, particularly those of high enthalpy, remain in thermo- 
dynamic and chemical equilibrium throughout. The purpose of this report is to extend 
the body of charts for carbon-dioxide flow to include those which can be used to esti- 
mate possible nonequilibrium effects in hypervelocity nozzles. It is difficult to obtain 
exact solutions of the nonequilibrium nozzle flow, and often the treatment of the atten- 
dant complexities is not justified in light of the results desired. In this study, esti- 
mates of nonequilibrium effects are based on the assumption of an instant freeze. 

This technique is described in more detail in reference 5. In this method it is 
assumed that the flow leaves the reservoir and travels in equilibrium to a certain 
point in the nozzle, beyond which the chemical and thermodynamic rate processes can 
no longer maintain equilibrium. At this nozzle position, termed the freeze point, it is 
assumed that the species concentrations remain constant (frozen) for the remainder of 



the expansion. The flow properties are then calculated with the molecular vibrational 
energy either remaining constant (frozen) at the value corresponding to the freeze point 
or continuing in equilibrium with translational and rotational temperatures. The prop- 
erties of flows in equilibrium are presented for comparison. Results of a similar 
study for air are charted in reference 6 . 

Results are presented for specified stagnation pressures of 1 to 1, 000 atm, stag- 
nation enthalpies of 2. 5 to 55 kj/g (1, 075 to 23,650 Btu/lbm), and assumed freeze 
points corresponding to Mach numbers of 1 through 6. The quantities presented are 
temperature, pressure, density, velocity, dynamic pressure, Mach number, Reynolds 
number, molecular weight ratio, and species concentrations. Temperature, pressure 
ratio, and density ratio across normal shock waves in a nozzle are also included. 


SYMBOLS 


a sound speed 

A area 

Cp specific heat at constant pressure 

c v specific heat at constant volume 

h enthalpy 

l characteristic length 

m molecular weight of mixture 

m G molecular weight of undissociated carbon dioxide 

M Mach number 

p pressure 

q dynamic pressure 

R universal gas constant 

Re Reynolds number 

T temperature 

u velocity 

x mole fraction, moles of species per mole of gas 

y mole fraction, moles of species per mole of undissociated gas at STP 


2 


molecular weight ratio, 


Z , 

7 


Po 

M- 

( ^s 
()* 


isentropic exponent, 


m o 
m 

5 Jn p\ 
3 In p/ 

^ /s 


density 

reference density, 1.979 kg/m 3 
coefficient of viscosity 
differentiation at constant entropy 
sonic point 


Subscripts 

f frozen condition 

0 reference condition 

t total condition 

1 condition in front of a normal shock wave 

2 condition behind a normal shock wave 


COMPUTATION PROCEDURE 


Nozzle-flow properties were calculated for three types of flow. In one type, 
the flow was assumed to be in chemical and thermodynamic equilibrium throughout. 

In the other two types, the flow was in equilibrium from the reservoir to points in the 
nozzle at which various assumed Mach numbers occurred. The chemical species con- 
centrations were evaluated at these points, and it was assumed that they remained 
frozen in these proportions for the remaining length of the nozzle. The two types of 
flow involving frozen chemistry differed in the manner in which the flow was expanded 
beyond the freeze point. In one case, the energy in molecular vibrations was held 
constant at the value consistent with the equilibrium conditions at the freeze point, and 
in the other, the vibrations were assumed to be in equilibrium with translational and 
rotational temperatures. 

The computational techniques are similar to those used in references 3 and 6, 
but the key steps are repeated here for completeness. All calculations were made 
with an electronic digital computer. 
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Equilibrium Flow 


A convenient method for calculating equilibrium flow quantities is to use the 
static pressure as the independent variable. It is desirable from the physical stand- 
point, however, to present results as functions of area ratio. The equilibrium ther- 
modynamic properties of reference 2 were used to determine the flow properties in 
the following manner. 

The equation of state, 


P m o 

ZRT 


and the energy equation, 


u = J2(ht - h) 


( 1 ) 


(2) 


form the basis for the computation. The thermodynamic properties of temperature, 
molecular weight ratio, and enthalpy were obtained from the results of reference 2. 

It was first assumed that there was an isentropic expansion from the reservoir pres- 
sure to a number of assumed values of static pressures, p. The product pu was then 
calculated for each value of the static pressure, and the maximum of this quantity was 
taken to be the mass flow per unit area at the throat. This value was then used in the 
continuity equation to calculate the area ratio, A/ A*: 


Mach number was obtained from 


A _ p*u* 
A* pu 


M 


u _ u 

a 7 ZRT/ m Q 


where 


a = speed of sound 



7 = isentropic exponent 



(3) 

(4) 


(5) 


( 6 ) 


The following equation was used to 
and unit length 


calculate the Reynolds number per unit pressure 

Re = UP_ 
li>t UPt 


0 ) 
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Values for the viscosity, p, of carbon dioxide were not available for most of the 
ranges of temperature and pressure considered. Therefore the viscosity was assumed 
to be the same as for air. Furthermore, the effects of dissociation and ionization on 
the viscosity were disregarded, so that viscosity could be calculated from Sutherland's 
formula for air (ref. 7) 

H = 1.09x1 0~ e T l/S + , &=§ (8) 

where T is in °K. Because of these approximations, emphasis should not be placed 
on quantitative results for Reynolds number but rather on the qualitative manner in 
which it varies with nozzle area ratio. 


The properties presented for the normal shock wave correspond to those that 
would occur if a normal shock wave existed at various values of area ratio. They were 
calculated in the following manner, with equilibrium flow assumed through the shock: 

1. The density ratio, p 2 /Pi , was assumed, and the pressure behind the shock 
was calculated from the combination of the momentum and continuity equations: 


p 2 = Pl + Pl u x 



( 9 ) 


2. The static enthalpy, h 2 , was calculated from the energy equation in the 
following form, 


h 2 = h-L + 


Ul 



( 10 ) 


3. Molecular weight ratio, Z 2 , and static temperature, T 2 , were then found 
from the results of reference 2 for the state properties p 2 and h 2 . 

4. Density, P 2 , was calculated from equation (1) and a new density ratio 
formed. 

5. The iteration was continued until two successive values of density ratio 
differed by a small amount (10 -4 ). 

The stagnation conditions behind the shock were calculated with the assumptions of an 
adiabatic process through the shock and an isentropic compression from the static to 
stagnation conditions. 

The species concentrations existing in the flow at various stations in the nozzle 
were obtained from reference 2 and were included so that frozen flows could be calcu- 
lated, as explained in the following section. 


Frozen Flow 

The flow was assumed to expand in equilibrium from the reservoir to a spec- 
ified location in the nozzle at which the chemical species concentrations were assumed 
to freeze instantaneously. This location, called the freeze point, was specified for a 
given equilibrium value of Mach number, denoted Mf, at that point. The subscript f 
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* 

identifies the equilibrium properties at the freeze point. The expansion beyond the 
freeze point was then assumed to occur with the energy in molecular vibrations either 
held constant at the value corresponding to the freeze point or continuing in equilibrium 
with translational temperatures. The two types of frozen flow are discussed in more 
detail in the following sections. 

Frozen chemistry and frozen molecular vibrations . - These calculations were 
made by first determining the specific heats, Cp and c v , and their ratio, 7, at the 
chosen freeze point. 


C P 


c v 


7 


YiCpi 

J 

(11) 

1 

yi c v, 

j 1 

(12) 

/ °v 

(13) 


The quantity, y i( denotes the concentration of the ith species in moles per mole of 
undissociated carbon-dioxide gas at STP (1 atm, 273° K). The species concentrations 
were determined from reference 2 for the equilibrium thermodynamic properties at 
the freeze point. The value for 7 in equation (13) is not the same as that for equilib- 
rium flow; consequently, the value of Mach number is discontinuous at the freeze 
point. 


With static pressure again the independent variable, the following equations 
were used. 





_p_ = /V\ /iYV 

p 0 W \p t y v p t^ 


h = h f + Cp(T - T f ) 


u = J2 [(h t -h f ) - c p (T - T f )] 

A = P*u* 

A* ~ PU 


a nT y ZRT/m 0 


(14) 


(15) 

(16) 
(17) 

(3) 

(4) 
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V(r-i) 


(18) 


(7 + 1)M‘ 




7 + 1 

-27M 2 - (7 - 1)- 


Equation (18) was used to determine the total pressure behind a normal shock 
wave with the assumption that the flow was frozen through the shock; whereas, in the 
equilibrium section this total pressure was found with the assumption of equilibrium 
flow through the shock. 

Frozen chemical reactions and equilibrium molecular vibrations. - The equa- 
tions of reference 1 are applicable to flows expanded from specified reservoir condi- 
tions with frozen chemistry and vibrational equilibrium. In the present study the flow 
with frozen chemistry and vibrational equilibrium is expanded from some specified 
freeze point (Mj ^ 1 ) downstream of the specified equilibrium reservoir condition. 
However, the equations of reference 1 may be used if appropriate frozen reservoir 
conditions can be found, such that the temperature and pressure at the freeze point 
are equal to the corresponding values for equilibrium flow. Since the chemistry is 
frozen, those quantities that depend on 7 will be discontinuous at the freeze point. 

The procedure for determining the appropriate stagnation conditions to use with 
the equations of reference 1 is as follows: 

1. At the freeze point, equilibrium values of temperature, species concentra- 
tions, and velocity were used with the equations of reference 1 to calculate 7 , speed 
of sound, and Mach number. 

2. To calculate the variation of stagnation temperature with static temperature 
the equations of reference 1 were used with the species concentrations and Mach num- 
ber determined in step ( 1 ). 

3. The appropriate stagnation temperature was determined from the results 
obtained in step ( 2 ) and the static temperature given in step ( 1 ). 

4. An expansion starting with the stagnation temperature of step (3) was cal- 
culated by the method of reference 1 ; however, only that portion beyond the freeze 
point was used. 

Note that the values of A* and total pressure associated with this frozen flow are 
different from those of the equilibrium flow upstream of the freeze point. The values 
of flow area and static pressure for the flow downstream of the freeze point were 
normalized with respective values of A* and p t for the equilibrium flow upstream of 
the freeze point. 
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RESULTS 


% 


All of the results are presented in graphical form in charts 1 through 18. An 
index to the charts is given in table I. Results are given for reservoir pressures of 
1, 100, and 1000 atm and reservoir enthalpies ranging from 2.5 to 55 kJ/g. Results 
obtained for flows with both chemistry and molecular vibrations frozen are presented 
for assumed freeze-point Mach numbers of 1, 2, 4, and 6. A limited number of cal- 
culations were made for flows with frozen chemistry and equilibrium vibrations, and 
results are presented for the quantities of Mach number, velocity, temperature, and 
pressure for a reservoir pressure of 100 atm and for reservoir enthalpies of 2.5, 20, 
and 55 kj/ g. Unit conversion factors are presented in table II to aid in converting 
results to other units. 
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TABLE I. - INDEX TO CHARTS 


Chart 

Ordinate 

Type of flow 

Page 

1 


T 

Equilibrium 

11 

2 

Temperature 

T 

Equilibrium and frozen 

12 

3 


T t 

i 2 

Equilibrium 

21 

4 


p/p t 

Equilibrium 

24 

5 

Pressure 

p/p t 

Equilibrium and frozen 

25 

6 


v\ 

Equilibrium and frozen 

31 

7 


(p/p 0 )/p t 

Equilibrium 

37 

8 

Density 

p/p 0 

Equilibrium and frozen 

38 

9 


P 2 /Pl 

Equilibrium 

44 

10 

Velocity 

u 

Equilibrium 

47 

11 


u 

Equilibrium and frozen 

48 

12 

Dynamic 

q/pt 

Equilibrium 

54 

13 

pressure 

q/p t 

Equilibrium and frozen 

55 

14 

..... . , 

Mach number 

1 

M 

Equilibrium 

57 

15 


M 

Equilibrium and frozen 

58 

16 

Reynolds number 

Re/ p t Z 

Equilibrium 

64 

17 

Molecular weight 
fraction 

Z 

Equilibrium 

67 

18 

Mole fraction 

X 

Equilibrium 

70 




TABLE II. - CONVERSION FACTORS 


Quantity 

Temperature 

Pressure 

Density 

Velocity 

Enthalpy 


Multiply 

By 

To obtain 

°K 

1.800 

°R 

N/m 2 

1.450X10' 4 

lb f/ in . 2 

N/m 2 

9.871X10' 6 

atm 

kg/m 3 

1. 940x1 0" 3 

slugs/ft 3 

m/s 

3.281 

ft/ s 

kJ/g 

4. 300x1 0 2 

Btu/lbm 
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Chart 1. - Variation of temperature with area ratio; equilibrium flow. 
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Chart 2. - Variation of temperature with area ratio; equilibrium and frozen flows. 
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A /A* 

(b) Pj. = 1 atm; h^, = 15. 0, 20. 0, 25. 0 kj/g 
Chart 2 . - Continued. 
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A/A* 


(d) p^. = 100 atm; h^. = 2. 5, 5.0, 10. 0 kj/ g 
Chart 2. - Continued. 
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Chart 2. - Continued. 
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A /A* 

(g) P t = 1000 atm; h t = 2. 5, 5.0, 10. 0 kj/g 
Chart 2. - Continued. 


18 




ioV 


T 



I0 1 



a /a* 

(h) P t = 1000 atm; h t = 15. 0, 20, 0, 25. 0 kj/g 
Chart 2. - Continued. 
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(a) = 1 atm 

Chart 3. - Variation of total temperature behind a normal shock wave with area ratio; 

equilibrium flow. 


21 




(b) p^. = 100 atm 
Chart 3. - Continued. 
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Chart 4. - Variation of pressure ratio with area ratio; equilibrium flows. 
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(a) p = 1 atm; h t = 2. 5, 5. 0, 10. 0, 15. 0, 20. 0 kJ/g 
Chart 5. - Variation of pressure ratio with area ratio; equilibrium and frozen flows. 




Chart 5. - Continued. 
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Chart 5. - Continued. 
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Chart 5. - Continued. 
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Chart 5. - Continued. 
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Chart 5. - Concluded. 
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(a) Pj. = 1 atm; hj. = 2. 5, 5. 0, 10. 0, 15. 0, 20. 0 kj/ g 

Chart 6. - Variation of total pressure ratio across a normal shock wave with area ratio; equilibrium and frozen flows. 
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Chart 6. - Continued. 
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Chart 6. - Continued. 
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Chart 6. - Continued. 
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Chart 6. - Continued. 
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Chart 6. - Concluded. 
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Chart 7. - Variation of density with area ratio; equilibrium flow. 
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(a) Pj. = 1 atm; h t - 2. 5, o. 0, 10. 0, 15. 0, 20. 0 kj/g 
Chart 8. - Variation of density with area ratio; equilibrium and frozen flows. 
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Chart 8. - Continued. 
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Chart 8. - Continued. 
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Chart 8. - Continued, 
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Chart 8. - Continued. 
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Chart 8. - Concluded. 







A/A* 

(b) p t = 100 atm 


Chart 9. - Continued. 
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Chart 10. - Variation of velocity with area ratio; equilibrium flow. 
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(f) p = 1000 atm; ht = 25. 0, 35. 0, 45. 0, 55. 0 kj/g 
Chart 11. - Concluded. 
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Chart 12. - Variation of dynamic pressure with area ratio; equilibrium flow. 
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A /A 

(a) h t = 2. 5, 5.0, 10.0, 15.0, 20. 0 kj/g 

Chart 13. - Variation of dynamic pressure with area ratio; equilibrium and frozen flows. 
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Chart 13. - Concluded. 



Chart 14. - Variation of Mach number with area ratio; equilibrium flow. 
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Chart 15 
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(b) p t = 1 atm; h t = 25. 0, 35. 0, 45. 0, 55. 0 kj/g 
Chart 15. - Continued. 
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(c) p - 100 atm; h t = 2.5, 5.0, 10. 0, 15. 0, 20. 0 kj/g 
Chart 15. - Continued. 
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(d) p = 100 atm; hj. = 25. 0, 35. 0, 45. 0, 55. 0 kj/ g 
Chart 15. - Continued. 
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(f) P t - 1000 atm; h t = 25. 0, 35. 0, 45. 0, 55. 0 kj/g 
Chart 15. - Concluded. 
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(a) Pj. - l atm 

Chart 16. - Variation of Reynolds number per unit length per unit total pressure with 

area ratio; equilibrium flow. 
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A /A* 

(b) p t = 100 atm 
Chart 16. - Continued. 


65 







d\/9d 



I0 1 I0 2 I0 3 I0 4 I0 5 

A/A* 

(c) } ^ 1000 atm 

Chart 16. - Concluded. 
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(b) Pj. = 100 atm 
Chart 17. - Continued. 
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(c) p = 1000 atm 
Chart 17. - Concluded. 
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(a) 2. 5 kj/g 

Chart 18. - Variation of species concentrations with area ratio; equilibrium flow. 
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Chart 18. - Continued. 
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Chart 18. - Continued. 
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Chart 18. - Continued. 
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Chart 18 . - Continued. 
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Chart 18. - Continued,. 
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Chart 18. - Continued. 





3 

o 

LO 

LQ 


% 


78 


NASA-Langley, 1965 


A-2055 


Chart 18. - Concluded 











